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Abstract 

The NA60 experiment at the CERN SPS has studied low-mass muon pairs in 158 
AGeV In-In collisions. A strong excess of pairs is observed above the yield expected 
from neutral meson decays. After subtraction of the decay sources, the shape of 
the resulting mass spectrum is largely consistent with a dominant contribution 
from 7r + 7r~ — » p — > annihilation. The associated p spectral function exhibits 

considerable broadening, but essentially no shift in mass. The acceptance-corrected 
Pt spectra have a shape atypical for radial flow. They also significantly depend 
on mass, pointing to different sources in different mass regions. Both mass and px 
spectra are compared to recent theoretical predictions. 
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1 Introduction 



Thermal dilepton production in the low-mass region is largely mediated by the 
light vector mesons p, u and 0. Among these, the p(770 MeV/c 2 ) is the most 
important, due to its strong coupling to the 7T7T channel and its short lifetime 
of only 1.3 fm/c, much shorter than the lifetime of the fireball. Theoretical 
predictions exist for changes both of the mass and the width which are tied 
directly [1,2] or indirectly [3] to chiral symmetry restoration. While broadening 
of the p is generally expected, changes of the mass are controversial. For masses 
>0.9 GeV/c 2 , other hadronic processes set in, and partonic processes start to 
play a larger role as well. The first NA60 results focused on the space-time 
averaged spectral function of the p [4]; more details were added recently [5]. 
The present paper therefore concentrates on new developments: an analysis 
of the shape of the excess mass spectra [5], and in particular first results on 
acceptance-corrected pt and my spectra for different mass windows. Mass and 
Pt spectra are also compared to the latest theoretical developments [6,7,8]. 



2 Experimental results 

Details of the NA60 apparatus are contained in [9,10], while the different anal- 
ysis steps (including the critical assessment of the combinatorial background 
from 7r and K decays through event mixing) are described in [11]. The results 
reported here were obtained from the analysis of data taken in 2003 for In-In 
at 158 AGeV. The left part of Fig. 1 shows the opposite-sign, background 
and signal dimuon mass spectra, integrated over all collision centralities. Af- 
ter subtracting the combinatorial background and the signal fake matches, the 
resulting net spectrum contains about 360 000 muon pairs in the mass range 
0-2 GeV/c 2 , roughly 50% of the total available statistics. The associated av- 
erage charged-particle multiplicity density measured by the vertex tracker is 
dN c h/dr] =120. The vector mesons u and <fi are completely resolved; the mass 
resolution at the u is 20 MeV/c 2 . A major part of the subsequent analysis is 
done in four classes of collision centrality (defined through the charged-particle 
multiplicity density): peripheral (4—30), semiperipheral (30—110), semicentral 
(110-170) and central (170-240). 

The peripheral data can essentially be described by the expected electromag- 
netic decays of the neutral mesons, i.e. the 2-body decays of the rj, p, uj and <f> 
resonances and the Dalitz decays of the 77, 77' and u [4,5]. This is not true in 
the more central bins, due to the existence of a strong excess. To isolate this 
excess with a priori unknown characteristics without any fits, we have used 
a novel procedure as shown in Fig. 1 (right), made possible by the high data 
quality. For each centrality bin, the cocktail of the decay sources is subtracted 
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Fig. 1. Left: Mass spectra of the opposite-sign dimuons (upper histogram), combina- 
torial background (dashed), signal fake matches (dashed-dotted), and resulting sig- 
nal (histogram with error bars). Rightdsolation of an excess above the hadron decay 
cocktail (see text). Total data (open circles), individual cocktail sources (solid), dif- 
ference data (thick triangles), sum of cocktail sources and difference data (dashed). 

from the data by using solely local criteria [4,5]. The p, on the contrary, is 
not subtracted. Examples for the resulting difference spectra are contained in 
Fig. 1 (right) for all centralities and in Fig. 6 for the semicentral bin. The 
qualitative features of the spectra are striking: a peaked structure is always 
seen, residing on a broad continuum with a yield strongly increasing with cen- 
trality, but remaining essentially centered around the position of the nominal 
p pole. Further details on the subtraction procedure and the results, including 
a discussion of the systematic errors (which reach up to 25% in the continuum 
region), can be found in [4,5]. 

A more quantitative analysis of the shape of the excess mass spectra vs. cen- 
trality has recently been performed, using a finer subdivision of the data into 
12 centrality bins. The data (Fig. 1, right) were subdivided into 3 windows, one 
central (C) in the mass range 0.64<M<0.84 GeV/c 2 , and two adjacent ones 
on the left (L) and the right (R) with equal widths. The spectral shape can 
then roughly be characterized by two yields: C-(L+R)/2 for the p-like peak, 
and 3(L+R)/2 for the underlying continuum. The ratios peak/p, continuum/p 
and peak/ continuum (where p stands for the cocktail p) are plotted in the left 
part of Fig. 2, compared to the RMS of the total mass interval 0.44<M<1.04 
GeV/c 2 on the right. Three essential features should be noted. First, the ratio 
peak/p decreases from the most peripheral to the most central bin by nearly a 
factor of 2, ruling out the naive view that the shape can simply be explained 
by the cocktail p residing on a broad continuum, independent of centrality. 
Second, the sum of the ratios (peak+continuum)/p is the total enhancement 
factor relative to the cocktail p; it reaches about 5.5 in the most central bin. 
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Fig. 2. Left: Yield ratios continuum/p, peak/p, and peak/continuum (see text). 
The errors shown are purely statistical; the point-to-point systematic errors are 
small compared to them. Right: RMS of the excess mass spectra in the window 
0.44<M<1.04 GeV/c 2 . The subtraction of charm follows the results from [11], with 
the exception of the most peripheral bin (100% charm at M=l.l GeV/c 2 ). 

Third, the centrality dependences are largely consistent in the two plots. In 
particular, the flat part up to dN c / l /dy=100 is followed by more rapid changes 
in all variables, statistically highly significant for the ratio continuum/p and 
the RMS. 



Dependences on pr were so far shown only as mass spectra associated with 
three different windows, without acceptance correction [5]. We have now 
obtained acceptance-corrected px spectra associated with three different mass 
windows. The NA60 acceptance relative to An as a function of pt is shown in 
Fig. 3 (left). In principle, the acceptance correction requires a 3-dimensional 




0.5 1 1.5 2 1-5 2 2-5 3 3.5 4 4.5 

p T (GeV/c) y 



Fig. 3. Left: Acceptance relative to 4-7T for different mass windows. Right: Rapid- 
ity distribution of the excess data for the mass window 0.6<M<0.9 GeV/c 2 and 
for three selected px bins. The measured data (full markers) are reflected around 
midrapidity (open markers). 
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grid in M — px — y space. However, the low populated bins in the phase 
space corners can introduce large errors once the correction is applied. To 
overcome this problem, the correction is performed in 2-dimensional M — p? 
space, using the measured rapidity distribution as an input. The latter was 
determined with an acceptance correction found, in an iterative way, from 
Monte Carlo simulations matched to the data in mass and pt- On the basis 
of this rapidity distribution, 0.1 GeV/c 2 bins in mass and 0.2 GeV/c bins 
in px were used to determine the remaining 2-dimensional correction. Once 
corrected, the results were integrated over the three extended mass windows 
0.4<M<0.6, 0.6<M<0.9 and 1.0<M<1.4. GeV/c 2 . The rapidity distribution 
of the central mass window is shown in Fig. 3 (right) for three different pt 
windows. It is interesting to note that it resembles closely the distribution of 
inclusive pion production, as measured by NA49 for Pb-Pb. 
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Fig. 4. Acceptance-corrected pt spectra for three mass windows and for three cen- 
trality bins. For discussion of errors see text. 
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Fig. 5. Left: Acceptance-corrected mj- spectra for three mass windows, summed 
over centralities (excluding the peripheral bin). Right: Inverse slope parameter T e ff 
from differential fits of the tut spectra (see text). 

The results for the acceptance-corrected px spectra are summarized in Fig. 4. 
The errors shown are purely statistical. Systematic errors arise from the ac- 
ceptance corrections including the rapidity distribution used, the subtraction 
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of the cocktail, and the subtraction of the combinatorial background plus fake 
matches. For pt<0.5 GeV/c, the combinatorial background contributes most, 
ranging from 10 to 25% for semiperipheral up to central. For p^>l GeV/c, the 
statistical errors dominate. The data show a significant dependence on mass, 
but hardly on centrality. To bear out the differences in mass more clearly, the 
data were summed over the three more central bins and plotted vs. trans- 
verse mass rriT in place of pt, as shown in Fig. 5 (left). The inverse slope 
parameter T e ff as determined from differential fits of the spectra with 
exp(— m,T/T e ff), using a sliding window in p T , is plotted on the right. The 
main features of the results are remarkable and somewhat unexpected. At 
very low all spectra steepen rather than flatten, as expected from radial 
flow, equivalent to very small values of T e ff (the resonance, placed just in 
between the upper two mass windows, does flatten as expected). Moreover, 
depending on the fit region, T e ff covers an unusually large dynamic range. 
Finally, the largest masses have the steepest m T spectrum, i.e. the smallest 
value of T e ff everywhere, again contrary to radial flow and to what is usually 
observed for hadrons. All this suggests that different mass regions are coupled 
to basically different emission sources. 
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Fig. 6. Comparison of excess mass spectrum, with no pt cut, to model predictions 
made for In-ln at dN c h j '(^=140 (semicentral bin) [6,7,8]. Data errors are purely 
statistical. The open data points show the difference spectrum resulting from a 
decrease of the i] yield by 10% (roughly the high-pr limit of the measured r\ yield). 
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3 Comparison to theoretical models 



The qualitative features of the excess mass spectra shown above are consistent 
with an interpretation as direct thermal radiation from the fireball, dominated 
by 7r + 7r~ — >■ p — > p + p~ annihilation. As explained in detail in [4,5], the ex- 
cess mass spectra, with no acceptance correction and no px selection, can 
be approximately interpreted as the spectral function of the intermediate p, 
averaged over momenta and the complete space-time history of the fireball. 
As already discussed in the first NA60 publications (see [4,5] and earlier ref. 
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Fig. 7. Comparison of acceptance-corrected mx spectra (left) and inverse slope 
parameters T e ff from differential fits (right) to model predictions made for In-ln 
at dN c /j/d?7=140 [6,7,8] (the data selection dN c / l /dry>30 corresponds to about the 
same value). Mass windows 0.6<M<0.9 (upper) and 1.0<M<1.4 GeV/c 2 (lower). 

cited ibd.), the moving-mass scenario a la Brown/Rho scaling [2] was found 
to be incompatible with the experimental excess mass spectra. Even consider- 
ing fireball variations within extremes, this has basically remained unchanged. 
Other scenarios, characterized essentially by a broadening of the p, have in 
the meantime been refined [6] or newly developed [7,8]. A summary of the 
comparison of these latest theoretical results to the excess mass spectrum for 
the semicentral bin is contained in Fig. 6. Hees/Rapp [6] have modified the 
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fireball part of the original predictions for In- In in the 2n region (based on [3]) 
and have now also added 47r processes, reflecting vector-axialvector mixing; 
the latter makes a strong contribution for M>1 GeV/c 2 . Renk/Ruppert [7] 
use the spectral function of Eletsky et al. including baryonic effects in the 2ir 
region, coupled to a strong partonic contribution which dominates the yield 
for M>1 GeV/c 2 . Dusling/Zahed [8] use a chiral virial approach including 
baryonic interactions, coupled again to some partonic contribution which is 
on a similar level as the hadronic yield for M>1 GeV/c 2 . All models achieve a 
reasonable description of the mass spectrum over the whole mass range, even 
in absolute terms. For M<0.5 GeV/c 2 , the decisive role played by baryonic 
interactions in [6] is seen in the upper right of Fig. 6 by switching them off. 
For M>1 GeV/c 2 , the agreement of all predictions with the data implies a low 
(if any) sensitivity to the dynamics here, equivalent to hadron-parton duality. 
The rriT spectra predicted by the models above are compared to the data in 
Fig. 7. In contrast to the mass spectra, none of the predictions is presently able 
to provide a satisfactory description of the my spectra over their full range. It 
seems that the rriT spectra are, in some sense, more sensitive to the dynam- 
ics than the mass spectra. A quantitative understanding of the spectra 
may thus help to discriminate between the different sources and provide fur- 
ther insight beyond hadron-parton duality (e.g. different radial flow effects for 
hadrons and partons). 

Summarizing, the previously measured excess mass spectra and the new accep- 
tance-corrected pr and uit spectra present rather strict boundary conditions 
to theoretical modeling. Beyond the p spectral function itself, this may also 
lead to a better understanding of the continuum part of the spectra. 
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